Searching for a Correlation Between Cosmic-Ray Sources Above 10^^ eV and 

Large-Scale Structure 



We study the anisotropy signature which is expected if the sources of ultra high energy, > 10 eV, 
cosmic-rays (UHECRs) are extragalactic and trace the large scale distribution of luminous matter. 
Using the PSCz galaxy catalog as a tracer of the large scale structure (LSS) we derive the expected all 
sky angular distribution of the UHECR intensity. We define a statistic, that measures the correlation 
between the predicted and observed UHECR arrival direction distributions, and show that it is 
more sensitive to the expected anisotropy signature than the power spectrum and the two point 
correlation function. The distribution of the correlation statistic is not sensitive to the unknown 
redshift evolution of UHECR source density and to the unknown strength and structure of inter- 
galactic magnetic fields. We show, using this statistic, that recently published > 5.7 x 10^'' eV Auger 
data are inconsistent with isotropy at ~ 98% CL, and consistent with a source distribution that 
traces LSS, with some preference to a source distribution that is biased with respect to the galaxy 
distribution. The anisotropy signature should be detectable also at lower energy, > 4 x lO'^^ eV. A 
few fold increase of the Auger exposure is likely to increase the significance to > 99% CL, but not 
to > 99.9% CL (unless the UHECR source density is comparable or larger than that of galaxies). 
In order to distinguish between different bias models, the systematic uncertainty in the absolute 
energy calibration of the experiments should be reduced to well below the current ~ 25%. 



The origin of UHECRs is still unknown: Sources have not been identified by cosmic-ray observations and all models 
of particle acceleration in known astrophysical objects are challenged by the extension of the spectrum to energies 
exceeding 10^° eV [H, The spectrum flattens (becomes harder) [3[ and there is evidence for a composition change 
from heavy to light nuclei [1, Q at ~ 10^^ eV. This suggests that while the cosmic-ray flux below ~ 10^^ eV is 
dominated by Galactic sources of heavy nuclei, it is dominated at higher energy by different sources of lighter nuclei. 
The isotropy of the UHECR arrival direction distribution suggests that these sources are extra- Galactic. 

If UHECRs are produced by a population of extra-Galactic astrophysical objects, that traces the distribution of 
luminous matter, the inhomogeneous distribution of luminous matter is expected to imprint a characteristic anisotropy 
signature on the UHECR arrival direction distribution. If UHECRs are indeed light nuclei, this anisotropy is expected 
to be significant above '--^ 5 x 10^^ eV, where the propagation distance of light cosmic-ray nuclei is limited by interaction 
with the background radiation field to a distance of the order of lO's of Mpc (d, 0, S|' see Q for a detailed review). 
This distance is comparable to that over which order unity variations in the density distribution of luminous matter are 
observed. Identification of the expected anisotropy signature would provide strong support to models where UHECRs 
are accelerated in known astrophysical objects, and would be inconsistent with most new physics, "top-down," models, 
in which UHECRs are produced by the decay of h eavy relic particles or topological defects [l| . 



The expected anisotropy signal was analyzed in [lO| . It was shown there that the signal should be detectable with 
high statistical significance once the number of detected UHECR events is increased beyond that available at that 
time, which was w 20 events above 4 x 10^^ eV, by a factor > 10. In this paper we revisit this subject, as the 
commissioning of the Auger detector may provide us within a few years with the required exposure [ll| . While the 
analysis of the current paper largely follows that of p^ . it is improved in several respects. First, for the derivation of 
the large-scale structure of luminous matter wc use a larger, as well as more complete and uniform, galaxy redshift 
survey, the PSCz catalog [l^ instead of the IRAS 1.2 Jy catalog The PSCz catalog contains 3 times more 

galaxies, including many more galaxies at distances of 150 — 300 Mpc. Second, we use an updated cosmological 
model, {O^ = 0.27, = 0.73, Hq = 75 km/s/Mpc} instead of {n^ = 1, J^a = 0, Hq = 100 km/s/Mpc}, and 
analyze the effects of a possible redshift evolution of the UHECR source density. Third, we study the sensitivity of the 
correlation statistic to systematic errors in the determination of UHECR energies. Finally, we compare our proposed 
method for identif ying the expected anisotropy signature to the more commonly used methods based on the power 
spectrum (e.g. [l^. Il5l. [TgI. [l7j) and on the two point correlation function (e.g. [H, [l^, [13, [2l[ and references therein) 
of the angular distribution of UHECR arrival directions. In order to facilitate the use of the proposed statistic for the 
detection of anisotropy in UHECR data, we have made available at http://www.weizmajin.ac.il/~waxman/criso 
a numerical description of the intensity maps (see ^ \V\ . 

An analysis of the expected anisotropy signal based on the PSCz catalog has recently been carried out in The 
most important effects taken into account in the present analysis, and neglected in ref. [2^ . are the dependence of 
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the anisotropy signal on the finite number density of sources and on the systematic errors in the determination of 
UHECR energies. As shown in ? IIIII and discussed in §|Vl both the finite number density of sources and the systematic 
uncertainties in energy determination have important impHcations to the predicted signal and its interpretation. 

We assume in our analysis that the UHERCs are protons. This assumption is motivated by two arguments. First, 
the observed spectrum of > 10^^ eV cosmic-rays is consistent with a cosmological distribution of proton accelerators 
producing (intrinsically) a power-law spectrum of high energy protons, d\ogn/d\ogE « —2, [23,[23| (see also fig. [T]). 
This intrinsic power-law spectrum is consistent with that expected in most models of particle acceleration ([3| and 
references therein). Second, the leading candidate extra-Galactic sources (GRBs and AGN, @| and references therein) 
are expected to accelerate primarily protons. 

The paper is organized as follows. In sec. [IT] we describe the formalism used to analyze the anisotropy signal and 
present maps of the predicted angular distribution of UHECR intensity. In sec. Illll we study the sensitivity of several 
statistics to the predicted anisotropy signature. In sec. IIVI we analyze Auger data, which was published during the 
preparation of this manuscript [25| . using our correlation statistic. Sec. fVl contains a brief summary of the results and 
a discussion of their implications. 

The following point should be made here regarding the isotropy analysis of ref. [2^. In ref. [l^l a correlation is 
found between the arrival direction distribution of > 5.7 x 10^^ eV cosmic-rays (detected by the Auger experiment) 
and between the angular distribution of low- luminosity AGNs included in the V-C AGN catalog [2^. According to the 
analysis of ref. [l^ , the probability that the detected correlation would arise by chance for an isotropic UHECR arrival 
direction distribution is ~ 1%. Since low luminosity AGNs trace the distribution of luminous matter, one may argue 
that the detected correlation provides the first evidence for a correlation of UHECR sources and LSS. Unfortunately, 
the V-C catalog is merely a compilation of AGN data available in the literature, and is therefore incomplete both in 
its sky coverage and in its luminosity coverage. It does not, therefore, provide a correct description of the local LSS, 
as clearly pointed out in the introduction of ref. [2^ : " The V-C catalog should not be used for any statistical analysis 
as it is not complete in any sense, except that it is, wc hope, a complete survey of the literature" . The interpretation 
of the correlation reported in ref. [l^l is thus unclear. 

II. METHOD 

We consider a model where the UHECR flux is produced by cosmological sources of protons tracing the large scale 
galaxy distribution. We assume that the sources are intrinsically identical and that the number density of sources 
is drawn from a Poisson distribution with an average given by 6[5]s(z), where s{z) is the average comoving number 
density of sources at redshift z and b is some bias functional of the local fractional galaxy over density, S = Sp/ p. The 
large scale structure galaxy density field is derived from the PSCz catalogue. We first derive in jj lll Al the 2D (angular) 
UHECR intensity map. The effect of inter-galactic magnetic fields on the UHECR arrival direction distribution, which 
is neglected in § III Al is discussed in detail in § IIIBI We show that deflections by the inter-galactic magnetic field 
do not affect significantly the UHECR intensity map on scales larger than a few degrees. In III CI we describe the 
Monte-Carlo method used to generate realizations of the cosmic-ray arrival direction distribution (The description 
given in III Al and III CI is brief, since it largely follows sections 2.1 and 2.2 of In ^ IIIDI we describe the statistics 

examined for the identification of the anisotropy signal. 

A. The cosmic-ray intensity map 

We estimate the large-scale galaxy density field by smoothing the galaxy distribution of the Point Source Catalog 
redshift survey (PSCz) [l^ with a Gaussian filter of variable dispersion given by the the greater of 6.4 h~^ Mpc and 
the mean galaxy separation [l^l (the number density of PSCz galaxies is weighted by the PSCz selection function [l^ 
before smoothing). This results in a smoothing length of 6.4 h~^ Mpc out to a distance of 50 Mpc, increasing 
to 10 Mpc at a distance of 100 Mpc, to 26 Mpc at a distance of 200 Mpc and to 50 Mpc at a 
distance of 300 h^^ Mpc. For distances greater than 300 Mpc we have assumed a homogeneous {Sp = 0) density 
field. 

The PSCz catalog covers about 84% of the sky, where most of the unobserved area is along the Galactic plane. 
This lack of coverage does not affect much the search for an extra-Galactic anisotropy signal. In such a search it 
is better to avoid using the arrival distribution of UHECR along the Galactic plane, since cosmic rays propagating 
through the Galactic plane may suffer strong deflections by the Galactic magnetic fleld. Nevertheless, in order to 
present complete intensity maps, we complete the density field in the unobserved regions following, e.g., [13, HI] by 
adding to the catalog "galaxies" with a number drawn from a Poisson distribution with average equal to the average 
density in the adjacent regions above and below the Galactic plane. 
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FIG. 1: The latest data on the highest energy cosmic rays. J is the differential intensity, J = dN/dEdAdtdfl. (a) The UHECR 
spectrum as published by the HiRes experiment [s^ and by the Auger experiment (preliminary, [ll[). (b) The spectrum after a 
shift of AE/E — +23% in the calibration of the absolute energy scale of the Auger experiment. The solid line is the spectrum 
that would be generated by a cosmological distribution of sources of protons, with intrinsic spectrum dlogn/dlogi? — —2 and 
redshift evolution following the star-formation rate, s{z) oc (1 + z)^. 



For the bias functional wc consider three models: an isotropic (I) model, in which the source distribution is uniform, 
b[6] = 1; an unbiased (UB) model where the source distribution traces the galaxy distribution with b[d] = 1 + 6; and 
a biased (B) model in which the UHECR source distribution is biased compared to the galaxy distribution with a 
threshold bias, b[S] = 1 + S for S > (Jmin and b[S] = otherwise. We use Smin ~ 0, which produces a source distribution 
which is more concentrated in the high density regions of the LSS (and hence also more concentrated towards the 
super-galactic plane) than IRAS galaxies are. 

As discussed in detail in [l^l, the local mean source density, s{z = 0), is constrained by the lack of repeaters 
to So = s(z = 0) > 10~^Mpc~'^. We will therefore consider in our analysis two values, sq = 10~^Mpc~'^ and 
So = lO^^Mpc^'', the latter corresponding to the number density of bright galaxies. The redshift evolution of the 
source density is also unknown. We therefore consider two extreme cases: no evolution, i.e. s(z) cx (1 + z)'^ , and 
fast evolution, s (x (1 -|- z)^, corresponding to the evolution of the star-formation rate [3l| and the AGN luminosity 
density [2^, isOl (the fastest evolution of known sources). We include in our analysis sources out to redshift 2 = 2 (the 
contribution of sources beyond that redshift to the flux above 10^^ eV is negligible). 

The intrinsic spectrum of protons produced by the sources is assumed to be dlogn/dlogE = —a with a w 2, 
extending to Smax = 10^^ eV. As mentioned in the introduction, this is the spectrum expected in models of particle 
acceleration in astrophysical sources. Furthermore, it was shown in [H, [23| that the UHECR spectrum produced 
by a cosmological distribution of sources of protons, with intrinsic spectrum dlogn/dlogE = —2, is consistent with 
the measured spectrum (as reported by the Yakutsk AGASA (s^l. Fly's Eye (s^l and HiRes experiments [11]). 
Figure [T] demonstrates that the preliminary spectrum of the Auger experiment is consistent with that of earlier 
experiments. The flux reported by the Auger experiment [llj is significantly lower than that reported by earlier 
experiments, e.g. HiRes, at energies ^ 10^^ eV. This discrepancy may be due to a systematic error in the absolute 
energy calibration of the experiments. In panel (b) of Fig. [T] we have adjusted the absolute energy calibration to 
bring the different measured fluxes into agreement at ^ 10^^ eV. This brings the measured flux into agreement at all 
energies > 10^^ eV (There is some discrepancy remaining at lower energies, probably indicating some errors in the 
calculated exposures, similar to the case for earlier experiments- see detailed discussion in [1^). The relative energy 
shift, \AE\/E = 23%, is well within the published systematic errors in energy scale. 

Figure [5] shows the contribution of sources at different distances to the differential flux of UHECRs at various 
energies (The contribution of sources within a given distance is lower than that presented in Fig. 1 of [loj since 
the current calculation assumes a rapid positive redshift evolution of source density). The figure indicates that at 
sufficiently large energy, > 50 EeV, the local, d < 100 Mpc, inhomogeneity of matter distribution should be imprinted 
on the angular distribution of UHECR intensity. 

Figures [3] and [4] present all sky maps of the UHECR intensity I{E, Cl), averaged over the realizations of the UHECR 



4 




FIG. 2: The fraction of differential cosmic-ray flux contributed by sources lying within a distance d from Earth. A homogeneous 
distribution of sources of protons, with intrinsic spectrum dlogn/dlogiS = — 2 and density evolving like the star formation 
rate, oc (1 + z)'^, was assumed. 

source distribution in the unbiased and biased models, for several UHECR energy thresholds E. I{E, (l) is given by 

Here, nQ{E) is the average number of protons of energy larger than E produced per unit volume and time at z = 0, 
and Eq{E,z) is the energy with which a proton should be produced at redshift z in order for it to be observed 
at z = 0, following energy loss due to interaction with the microwave background radiation, with energy E (We 
calculate Eq{E,z) using the continuous energy loss approximation, following [23}). The angular structure of the 
UHECR intensity map reflects the local large scale structure of the galaxy distribution, as can be seen by comparing 
figures [3] and m to figure O which presents the integrated galaxy density out to a distance of 75 Mpc. One may clearly 
identify the "Great Attractor", composed of the Hydra-Centaurus (300° < I < 360°, 0° < 6 < 45°) and Pavo-Indus 
(320° < / < 360°, -45° < 6 < 0°) super clusters, the Perseus-Pisces super cluster (130° < I < 160°, -30° <b< 30°), 
the NGC1600 group {I ~ 200°, b 30°), and the local voids along the Galactic plane (at 0° < / < 120°). 
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B. Magnetic fields 

UHECR may suffer significant deflections as they cross dense large scale structures, like galaxy clusters and large 
scale galaxy filaments, in which the energy density of the plasma is large enough to support strong magnetic fields. 
Let us consider first deflections by massive clusters of galaxies. As pointed out in [3^, observed magnetic fields are 
consistent with the assumption that the magnetic field energy density amounts to ~ 1% of the plasma thermal energy 
density (i.e. 1% of equipartition). For a typical temperature of a few keV and density of 10~'*cm~'^ over a scale of 
1 Mpc, this gives B ~ 0.1 /iC, consistent with observations [13, [H, [3^. The deflection of an UHECR passing through 
a cluster is approximately given by [4^ 

2\ {DXy/^- ^ /_L ^\ /^x 1/2 f_^\ ' 

9 J E/eB ■ VlMpclOkpcy VO.Ol/ \m^Oe\ J ' 




FIG. 3: All sky maps of the UHECR intensity I{E,Q), averaged over the realizations of the UHECR source distribution in 
the unbiased model (eq. [T]with b[S] = 1 + 5), for several UHECR energy thresholds E. The coordinates are Galactic and / 
is normalized to its all sky average, I = J dill (Q) / An . Model parameters (intrinsic source spectrum dlogn/dlogjS = —2 and 
redshift evolution cx (1 + z)^) are similar to those of Figure[2] The contours denote I/I — (0.7, 0.9, 1, 1.1, 1.3, 1.5), with dashed 
lines representing under-density. The solid green line denotes the super-galactic plane. The dashed-dotted blue line marks the 
boundary of Auger's coverage (corresponding to a zenith angle of 60°). 



Biased model 




FIG. 4: Same as figure [S] for the biased (B) model, b[5] = 1 + 5 for 5 > 0, b = otlierwise. Tiie contours denote 
(0.33,0.67,1,1.5,2,2.5,3). 
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FIG. 5: The integrated galaxy density out to a distance of 75 Mpc, normalized to the mean integrated density. The contours 
are logarithmic, ranging from 0.5 to 4 with three contours per density doubling. Dashed curves represent under-density. 



where es is the ratio of magnetic field to thermal plasma energy density {B ~ -v/STresnT), L is the propagation distance 
through the magnetized plasma, and A is the field correlation length, which is typically inferred from Faraday rotation 
measures to be of order 10 kpc [HI, S^]. Thus, a deflection of a few degrees is possible if the UHECR passes through 
a massive cluster. However, massive clusters are very rare objects, which contain a negligible fraction of the mass in 
the universe and cover only a negligible fraction of the sky. 

Passing through a large scale filament is more likely, but produces much smaller deflections: For a typical density of 
~ 10~^cm~'^ over a scale of 1 Mpc and temperature of ~ 0.1 keV [i^, the magnetic field is expected to be (assuming 
1% of equipartition) i? ~ 5 nG leading to a deflection 

1/2 



^^01° ^] (^] { ^ \ 

1 Mpc 10 kpc y VO.Ol/ \W^(N ) 



(3) 



Assuming that a fraction / ~ 0.1 of the volume is filled with filaments of diameter L ^ \ Mpc, the number of filaments 
crossed by an UHECR propagating a distance o? is ~ fd/L, and the total defiection is thus 



1 Mpc V 0.1 100 Mpc 10 kpc y VO.Ol/ V102"eV 



[^)' . (4) 



The results derived above are consistent with the recent semi-analytic analysis of [4J], who obtain 6* w 3° for 
E = 10^° eV and d = 100 Mpc, assuming = 0.1 and A = 100 kpc. The values adopted in Q for and A, 
and hence the deflection 6 derived in [ij, should be considered as upper limits: As noted above, cb = 0.1 is a value 
larger than typically inferred, and building an equipartition field with A = 100 kpc coherence scale is marginal since 
the turn over time of A = 100 kpc eddies in the expected LSS flow is comparable to the Hubble time. Our results 
are also consistent with, and provide a simple explanation of, the results of detailed numerical calculations. Dolag et 
al. [4^ have used numerical simulations of large scale structure formation, constrained to reproduce the large scale 
structure determined by the IRAS catalog out to ~ 100 Mpc, in order to study the expected deflection of UHECR by 
inter-galactic magnetic flelds. Assuming that the magnetic fleld energy density constitutes a signiflcant fraction of the 
(turbulent and thermal) energy density of the inter-galactic plasma, they conclude that the deflection of > 40 EeV 
UHECR does not exceed 2 degrees over 99% of the sky. 

The deflection of 4x 10^^ eV protons over d ^ 100 Mpc propagation distance is therefore expected to be smaller than 
or comparable to a few degrees. Note, that although the deflection may be somewhat larger for longer propagation 
distance, the largest propagation distance of > 4 x 10"'^^ eV protons is ^ 350 Mpc (see fig.H]), the random deflections 
of particles originating from sources at distances > 100 Mpc are not expected to affect significantly the intensity 
distribution: The anisotropy of the angular UHECR intensity distribution is dominated by structures at distances 
d < 100 Mpc, while sources at larger distances produce a roughly isotropic "background". We therefore expect 
deflections not to modify the intensity distribution of > 4 x 10^^ eV on scales larger than a few degrees. 

It should be pointed out here that large distortions of the UHECR intensity map could have been generated by 
inter-galactic magnetic fields had the Milky Way galaxy been embedded in a high density large scale filament or 
cluster, which may support a large magnetic field. However, both observations and constrained simulations of our 
local super-cluster, e.g. (43| . show that the Galaxy does not lie within or near such a structure. The discrepancy 
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between the results of [45| and those of [46[, who predict much larger deflections based on unconstrained large scale 
structure simulations, may largely be due to the presence of highly magnetized large scale structures near, or around, 
the Galaxy, which are inconsistent with the observed nearby galaxy distribution (see [i^ for detailed discussion). 



Monte-Carlo realizations 



Following |10 | we denote by Si{E,Cl) the number of sources per unit solid angle in the direction Cl, that produce 
i events in the detector above energy E. It was shown in [lo| that Si and Sj are statistically independent for j ^ i, 
and that 5*^ is Poisson distributed with an average 



S, {E,n) = J dzc 



SizM) P,{E,z). (5) 



Here, is the luminosity distance and Pi{E, z) is the probability that a source at redshift z produces i events above 
energy E. Pi is given by a Poisson distribution, 

P,(i?,.) = ffii^e-^(^-), (6) 

where N{E^ z) is the average number of events above energy E produced by a source at redshift z. N{E, z) is given by 
the product of the effective detector area, A, the observation time, T, and the flux F{E, z) of particles above energy 
E produced by a source at redshift z, N{E, z) = F{E, z)AT. F{E, z) is related to the rate N{E) with which particles 
of energy larger than E are produced by a single source by [To| 

Here, Eo{E, z) is the energy with which a proton should be produced at redshift z in order for it to be observed at 
z = 0, following energy loss due to interaction with the microwave background radiation, with energy E. Expressing 
N{E) as a ratio between ho{E), the average number of protons of energy larger than E produced per unit volume 
and time at z = 0, and sq, we finally have [lO| 

So 4:TrdL(zy 

Note, that the effective area and observing time of a particular experiment, AT, may in general depend on Cl, in 
which case N and Pi become Cl dependent as well, N{E, z, Cl) = F{E, z)AT{Cl). 

Using eqs. ([5]), ([6|) and ^ we generate a realization of the cosmic-ray arrival direction distribution by drawing the 
number of sources producing i events above energy E in an angular region dQ around direction from a Possion 
distribution with average Si{E, fl)dfl. The production rate ho is normalized based on the observed UHECR flux. The 
normalization is uncertain due to the uncertainty in the absolute energy calibration of the UHECR experiments ([2^. 
see also fig. [Ij. Adopting, for example, the absolute energy scale of the HiRes experiment, the average number of 
events detected above 40 EeV is 300 per year for the Auger exposure of 7000 km^ sr, and about 50 events per year for 
the planned Telescope Array exposure of 1200 km^ sr [49|. These numbers are reduced by a factor of approximately 
2 if the preliminary absolute energy scale of the Auger experiment is chosen (see fig. [T|), and increased by a factor of 
about 1.5 if the AGASA absolute energy scale is adopted (see fig. 1 of [23|). 

We take into account the statistical uncertainty in the energy determination of the experiment by modifying N{E, z): 
N{E,z) is replaced with J dE'[dN{E' , z)/dE']g{E' , E) where g{E',E) is the probability that an event of energy E' 
will be measured to have energy exceeding E. This procedure correctly reproduces the average observed number of 
events, N{E,z), but leads to some overestimate of the variance of the observed number of events. For a relative 
statistical energy uncertainty /S.E / E, we assume that the measured energy E follows a Gaussian distribution around 
the true energy E' , with a standard deviation of AE{E'). As we show in § IIIIl the statistical uncertainty in energy 
determination docs not affect the results significantly. 



D. Measures of anisotropy 



We consider several methods for the detection of the anisotropy of the cosmic-ray arrival distribution, using statistics 
based on the two point correlation function, on the angular power spectrum, and on a correlation between the observed 
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and predicted arrival direction distributions. The various statistics are described below. In § IIIII we compare the 
sensitivity of these statistics to the expected anisotropy signal by deriving, using the Monte Carlo simulations, their 
expected distributions for the different models (isotropic, unbiased, and biased UHECR source distribution). 

The (cumulative) two-point correlation function is defined as the number of pairs of events that have a distance 
smaller then an angular distance D (e.g. [i?}). 

N 

w{D) = J2Y.e(^-D.,,) (9) 

i j<i 

where Q is the step function, TV is the number of CR events, and Dij is the angular distance between CR events i 
and j. We use several statistics based on the two-point correlation function, all of the form 

y ur.^^-^ iW{D^, A+i) ~ WUD„ A+i)]' 

Xw[{L),^) = 1., , (10) 

(T^(A, A+i) 

where {Di} is a set of angular distances, Di < Di+i, W{Di, Di^i) = W{Di^i) — W{Di), Wiso{Di, Di^i) = 
Wiso{Di+i) — Wiso{Di), Wiso{D) is the average value of the correlation function obtained in the isotropic model, 
and a^{Di,Di+i) is the variance of Wiso{Di, Di+i) obtained in the isotropic model. 

The 2D angular UHECR intensity map may be decomposed into spherical harmonics, = y^,e af,nYp,m(tl). 

The observed UHECR arrival direction distribution may be used for estimating the coefficients agm by [i^ 

1 ^ 1 

ai,n = — Y,—Yirn{n'). (11) 

z— 1 

Here N is the number of discrete arrival directions {^2'}, Af = ^-nd Wi is the relative experimental exposure 

at arrival direction Cli (e.g. [l^). The angular power spectrum coefficients are given by 



Ci = — 

^ 21 



qry E (12) 



We use several statistics based on the two-point correlation function, all of the form 

XK(m) = s,^^i-§^, (13) 

where {£} is a set of £ values, Ciso,e is the average value of Ci obtained in the isotropic model, and af is the variance 
of Ci obtained in the isotropic model. 

Finally, we consider a statistic characterizing the correlation between the predicted and observed arrival distribu- 
tions, 

Xc,M - 2^ — — ■ ■ . (14) 

( .-, ^ ^ i.iso 

{«} 

Here {i} is a set of angular bins, iV^ is the number of events detected in bin i, Nij^o is the average number of events 
expected to be detected in an isotropic model, and Ni^M is the average number of events expected to be detected in the 
M (e.g. unbiased, biased) model. Throughout the paper we calculate Xq using 6° x 6° angular bins, in order to avoid 
sensitivity to the possible distortion of the arrival direction distribution by inter-galactic magnetic field deflections 
(see s ura . 



III. RESULTS 



Let us first compare the sensitivity of the various statistics to the expected anisotropy signal. For this comparison 
we adopt a local (z = 0) source density sq = 10~^Mpc~'^, a source density evolution following the star- formation 
rate, s oc (1 -I- z)^, and an intrinsic source spectrum d log n/d log i? — —2. We consider the angular arrival direction 
distribution of 100 cosmic-rays above 40 EeV, detected by an experiment with a statistical uncertainty of (single) 
event energy determination of AE/E = 20%. We consider a (hypothetical) detector with uniform exposure over the 
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FIG. 6: Properties of the distributions of Xc, W{Di, -Di+i) and Ce obtained from 10000 Monte-Carlo realizations of 300 events 
above 4 x 10^^ eV for the isotropic model (blue diamonds), unbiased model (green triangles, b[S] = 1 + 5) and biased model 
(red squares, b[S] = 1 + S for 5 > 0, b[5] = otherwise): (a) Probability distribution of Xc,ub', (b) Average and 1 a range of 
W{Di, -Di+i) — Wiso{Di, -Di+i) as a function of Di+i; (c) Average and 1 a range of Ce — Ciso,e as a function of £. We used here 
So = 10~*Mpc~^, s cx (1 + z)'^ , dlogn/dlog E = —2, AE/E = 0.2, and a uniform sky exposure. 



sky, and a detector with an angular ex pos ure pattern corresponding to that of the Auger experiment (as given in |14l |) 
or of the Telescope Array experiment [49[. For the current comparison we do not exclude the Galactic plane strip. 

Fig. [6] presents properties of the distributions of Xc,ub (cq. [14]), W{Di, Di+i) (eq. ^ and C; (eqs. [TT] & [T2l) . 
obtained from 10000 Monte-Carlo realizations of the isotropic, unbiased and biased models. We compare in table U 
the probability distributions of Xq^ub, Xq^b, Xy{{£}), and XwUDi}) for several sets of {£} and {Di} where the 
signal is strongest (/ < 4 and D < 60°). The table presents the probabilities P{Mi\M2) to rule out model Mi with 
confidence y, assuming the UHECR source distribution follows model M2: P{Mi\M2) is the probability to obtain, 
assuming model M2, a value of X {Xc, Xy, Xw) that is obtained in model Mi with probability smaller than 1 — y. 
Results are presented in the table for the detection of 100 cosmic-rays above 4 x 10^^ eV. Clearly, the Xc statistic 
provides the highest probability to detect the expected anisotropy signal. 

Next, we consider the dependence of the strength of the anisotropy signal on the cosmic-ray energy threshold and 
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95% CL 






99% CL 






P(I|UB) [%] 


P(I|B) [%] 


P(UB|B) [%] 


P(I|UB) 


%] P(I|B) [%] 


P(UB|B) [%] 


Xc,uB Auger 


23 


79 


42 


5 


43 


10 


Xc,B 


26 


73 


32 


7 


40 


6 


Xc,uB Full Sky 


22 


83 


48 


6 


51 


19 


Xc,B 


26 


78 


37 


7 


45 


9 


Xc.uB TA 


21 


82 


48 


5 


46 


11 


Xc,B 


25 


78 


37 


5 


37 


10 


Xwi{Di = 20,60}) Auger 


6 


9 


7 


2 


3 


1 


Xw{{D, = 10,20}) 


6 


11 


9 


2 


4 


2 


A = 0,10,20,30,40}) 


7 


12 


10 


2 


4 


2 


Xwi{D^ = 20,60}) Full sky 


6 


8 


7 


1 


1 


1 


Xwi{D, = 10,20}) 


4 


7 


8 


1 


1 


1 


Xw{{D, = 0,10,20,30,40}) 


5 


8 


7 


1 


1 


2 


XyHI} = {1}) Auger 


5 


6 


6 


1 


1 


1 


Xy{{£} = {2}) 


6 


8 


7 


1 


2 


1 


Xy(W = {1,2,3,4}) 


6 


7 


6 


1 


1 


1 


Xviie} = {1}) Full sky 


6 


10 


8 


1 


2 


2 


Xvm = {2}) 


6 


10 


8 


1 


2 


2 


Xy(W = {1,2,3,4}) 


7 


11 


7 


1 


2 


2 



TABLE I: Probabilities P{Mi\M2) to rule out model Mi at a certain confidence level (CL), assuming that the UHECR source 
distribution follows model M2, for a given statistic (Xc , Xy , Xw)- Numbers in each line are calculated from 10000 Monte Carlo 
realizations of 100 cosmic-rays above 40 EeV. Boxes marked "Full sky" refer to a uniform exposure over the sky, and boxes 
marked "Auger" ("TA") refer to an angular exposure pattern corresponding to that of the Auger (Telescope array) experiment. 
For this comparison we have used so = lO'^Mpc""^, s oc (1 + z)^, dlogn/dlog E = —2, and AE/E = 0.2, and included events 
along the Galactic disk. 

on the statistical energy uncertainty. Table [III presents the probabilities P{Mi\M2) for several energy thresholds, 
E > 20, 40, 60, 80 EcV, and parameters similar to those chosen for fig. [6l s cx (1 + z)^, d\ogn/d\ogE = —2, 
AE/E = 0.2 and angular exposure pattern corresponding to that of the Auger experiment. Probabilities are given 
for both So = lO^^'Mpc"'^ and sq = 10~^Mpc~^. The probabilities were calculated for this table excluding the 
Galactic strip \b\ < 12°, along which large deflections are expected by the Galactic magnetic field. Comparison of 
tables [IT] and [T] demonstrates that excluding the Galactic strip does not affect the signal significantly. The table 
demonstrates that the signal is stronger for lower energy. Although the contrast of the fluctuations in the cosmic-ray 
intensity is higher at high energy (see fig. [3]), due to the fact that the propagation distance is smaller, the signal 
becomes weaker at higher energies since the number of observed UHECR drops rapidly with energy. At £^ > 20 eV a 
significant contamination from Galactic sources may be present and the deflection of cosmic-ray particles may become 
significant. In order to avoid distortions of the UHECR intensity map by these effects, it is advisable to chose an 
energy threshold of > 40 EeV. The last two lines of each block of table [TTl demonstrate that statistical energy 
uncertainty of AE/E < 0.2 does not affect significantly the results. 

A comparison of the values of P{Mi\M2) given in table |TT] for sq = lO^^Mpc"^ and sq = lO'^Mpc"^ shows that 
the interpretation of the anisotropy signal (i.e. of the measured value of Xc) depends on the assumed source number 
density. While the average value of Xc obtained in realizations of the various models (I, UB, B) is independent of 
So, its distribution does depend on sq: For a lower source number density the number of sources contributing to the 
flux is lower, hence fluctuations are larger and the distribution of Xc values is wider. This is illustrated in panel c 
of fig. [71 Thus, in order to determine the probability of obtaining a certain value of Xc in a given model (I, UB, 
B), one must make an assumption regarding the value of so- As mentioned in § [TTl ^-nd discussed in detail in p^ . 
s(z = 0) may be constrained by the number of "repeaters" , the number of sources producing multiple events. The last 
column of table [TTl presents the probability for observing repeaters at various energies. Comparison of the numbers 
obtained for sq = 10~'*Mpc^'^ and sq — 10~^Mpc~'^ shows that the presence or absence of repeaters may provide 
additional constraints on the source density (beyond sq > lO^^Mpc"^) once the number of events observed above 
40 EeV exceeds several hundreds. 

The values of P{Mi\M2) that appear in table [TTl demonstrate, in accordance with the results of [l^l, that the 
detection of few hundred events above 4 x 10^^ eV is required in order to guarantee identification of the predicted 
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95% CL 








99% CL 






99.9% CL 






P{Mi\M2) [%] 


P(I|UB) 


P(I|B) 


P(UB|B) 


P(I|UB) 


P(I|B) 


P(UB|B) 


P(IIUB) 


P(IiB) 


P(UB|B) 


Pr 


E>20 EeV (402 events) 


26 


(31) 


67 (76) 


27 


(34) 


7(12) 


34 (51) 


7(14) 


<1 (3) 


1 (25) 


<1 (3) 


99 (16) 


E > 40 EeV 100 events 


22 


(26) 


67 (76) 


30 


(38) 


5 (10) 


30 (52) 


7(16) 


<1 (2) 


<1 (24) 


<1 (4) 


65 (2) 


E>60 EeV (31 events) 


20 


(25) 


60 (73) 


26 


(36) 


4(9) 


23 (48) 


6 (15) 


<1 (2) 


<1 (19) 


<1 (4) 


31 (1) 


E>80 EeV (10 events) 


15 


(19) 


38 (50) 


17 


(22) 


3(6) 


11 (25) 


3(7) 


<1 (1) 


<1 (8) 


<1 (1) 


9 (0.1) 


E>40 EeV AE/E = 0.1 


22 


(27) 


69 (80) 


33 


(41) 


5 (10) 


32 (56) 


9 (18) 


<1 (2) 


<1 (28) 


<1 (5) 


65 (2) 


E >40 EeV AE/E = 


23 


(27) 


70 (79) 


33 


(42) 


5 (10) 


33 (56) 


8 (18) 


<1 (2) 


<1 (28) 


<1 (5) 


64 (2) 


£ > 20 EeV (1205 events) 


45 


(62) 


94 (99) 


45 


(67) 


13 (36) 


70 (96) 


11 (41) 


<1 (12) 


3 (82) 


<1 (16) 


>99.9(90) 


E > 40 EeV 300 events 


39 


(54) 


94 (99) 


52 


(74) 


10 (29) 


70 (96) 


15 (49) 


<1 (10) 


3 (84) 


<1 (23) 


>99.9(37) 


£ > 60 EeV (94 events) 


31 


(50) 


87 (99) 


42 


(72) 


6 (25) 


44 (94) 


9 (45) 


<1 (7) 


2 (78) 


<1 (17) 


97 (15) 


E >80 EeV (31 events) 


22 


(36) 


63 (89) 


24 


(46) 


3 (16) 


14 (71) 


4(21) 


<1 (4) 


<1 (39) 


<1 (6) 


65 (4) 


E>40 EeV A£/S = 0.1 


40 


(55) 


96 (99.5) 


57 


(79) 


10 (29) 


74 (97) 


17 (54) 


<1 (10) 


5 (87) 


<1 (26) 


>99.9(35) 


E>40 EeV AE/i; = 


39 


(55) 


96 (99.5) 


59 


(80) 


10 (29) 


74 (97) 


19 (55) 


<1 (9) 


4 (87) 


<1 (26) 


>99.9(35) 


E >20 EeV (2410 events) 


62 


(86) 


99.4 (99.9) 


63 


(90) 


27 (66) 


94(>99.9) 


24 (72) 


3 (33) 


46 (99.4) 


2 (44) 


>99.9(>99.9) 


£ > 40 EeV 600 events 


52 


(78) 


99(>99.9) 


67 


(94) 


14 (54) 


88(>99.9) 


21 (81) 


<1 (26) 


5 (99.6) 


<1 (55) 


>99.9(88) 


E >60 EeV (188 events) 


41 


(74) 


95(>99.9) 


52 


(93) 


7 (49) 


61 (99.9) 


12 (78) 


<1 (21) 


2 (99.1) 


<1 (48) 


>99.9(60) 


E >80 EeV (62 events) 


27 


(57) 


76 (99.2) 


29 


(69) 


3 (30) 


21 (95) 


5 (42) 


<1 (9) 


1 (79) 


<1 (16) 


97 (24) 


E>40 EcV AB/S = 0.1 


53 


(80) 


99.2(>99.9) 


73 


(96) 


13 (57) 


90(>99.9) 


23 (86) 


<1 (27) 


5 (99.7) 


<1 (62) 


>99.9(87) 


E >40 EeV A£/£; = 


54 


(80) 


99.3(>99.9) 


74 


(97) 


14 (56) 


91(>99.9) 


24 (87) 


<1 (27) 


5 (99.8) 


<1 (62) 


>99.9(87) 



TABLE IL The energy threshold and energy resolution dependence of P{Mi\M2) for the Xc,ub statistic. Model parameters 
used: s oc (1 + z)^ , d\ogn/d\og E — —2, and AE/E = 0.2. Values are given for sq = 10~*Mpc~^ and in parentheses for 
So — 10~^Mpc~"^. An angular exposure pattern corresponding to that of the Auger experiment was assumed, and the Galactic 
strip |fe| < 12° was excluded. Numbers in each line are calculated from 100,000 Monte Carlo realizations. The top block 
of lines gives the results for 100 cosmic-rays above 40 EeV, the middle block for 300 events, and the bottom block for 600 
events (numbers in parentheses in the first column give the corresponding numbers of events above different energy thresholds) . 
The last two lines in each part of the table illustrate the effect of statistical energy uncertainty. The last column gives the 
probability for the detection of "repeaters", sources producing multiple events, for so — 10~*Mpc~'^ and in parentheses for 
So = 10-^Mpc-^ 

anisotropy signature, and discrimination between different bias models, with high statistical significance. 

Finally, we consider the sensitivity of the Xc.ub distribution to variations in dlogn/dlogE, to modifications of the 
redshift evolution of the source number density, and to uncertainties in the calibration of the absolute energy scale 
of the experiments. The various panels of fig. [7] compare the distributions of Xc,ub obtained for sq = 10~'*Mpc~'^, 
s oc (1 + zY, d\ogn/d\ogE = —2, AE/E = 0.2 and £' > 40 EeV with those obtained with modified parameters 
(for all calculations we have used an angular exposure pattern corresponding to that of the Auger experiment and 
excluded the Galactic strip |6| < 12°). Panels (a) and (b) demonstrate that the distribution of Xc.ub is not sensitive 
to variations in dlogn/dlogE and in the redshift evolution of the source number density. Panel (c) addresses the 
sensitivity to variations of sq, which was already discussed above. 

Panel (d) addresses the sensitivity to uncertainties in the calibration of the absolute energy scale, ft compares the 
distributions of Xc,ub which are obtained by calculating Xc,ub (eq. I14p with the model average number of events 
expected for E > 40 EeV, Ni^usiE > 40 EeV), and the number of events obtained in the Monte-Carlo simulations 
for different energy thresholds, Ni{E > E') with E' = 30 EeV, 40 EeV, and 50 EeV (properly normaHzed). The 
curve obtained for E' = 30 EeV (50 EeV) simulates the distribution of Xc values that are expected to be obtained 
for a systematic overestimate (underestimate) of event energies by the experiment. The figure demonstrates that 
the distribution of Xc.ub values, which is expected to be measured from the distribution of UHECR events with 
energies inferred hy the experiment to satisfy E' > 40 EeV, is sensitive to uncertainties of 25% in the experimental 
calibration of the absolute energy scale: Since the anisotropy is stronger at higher energy, a systematic underestimate 
(overestimate) of event energies leads to a shift of the Xc distribution towards larger (smaller) values. Since the 
anisotropy is also stronger (at a given energy) for stronger bias, this implies that an underestimate (overestimate) of 
the absolute energy scale will lead to an overestimate (underestimate) of the bias of the UHECR source distribution. 
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(a) (b) 




FIG. 7: The sensitivity of the distribution of Xc,ub to variations in model parameters and to systematic energy uncertainty. 
All curves obtained from 10000 Monte-Carlo realizations with 300 events above 40 EeV. In each panel, solid lines correspond 
to default model parameters, so — 10~''Mpc~'^, s oc (1 + z)"^, dlogn/dlog E = —2, and blue/green/red color corresponds 
to the isotropic/unbiased/biased model. Panels (a-c) explore the following variations: (a) Dashed lines correspond to a = 
dlog n/dlog -E = —2.2; (b) Dashed lines correspond to no source evolution, s cx (1 + (c) Dashed lines correspond to 
So = lO^^Mpc^'*. In panel (d) the sensitivity to systematic energy uncertainty is addressed by calculating Xc,ub fea.ll4p with 
Ni^UB^E > 40 EeV) and Ni obtained from the simulations for different energy thresholds (and properly normalized): 30 EeV 
(dashed line), 40 EeV (sohd line), and 50 EeV (thin-solid line). 

IV. ANALYSIS OF > 5.7 x lO^'-* eV AUGER DATA 

Figure [5] shows the positions of the 27 Auger events with energy exceeding 5.7 x 10^^ eV, overlaid on the intensity 
map obtained in the biased model (see figure S]) . Examination by eye indicates an enhancement of the UHECR flux 
associated with the Great Attractor. The significance of this apparent correlation is quantified in table IIIll The table 
presents the value of Xq.ub obtained for the Auger data, the average values oi Xc,ub that are expected to be obtained 
in the various models (isotropic, unbiased, biased), and the probability that the value oi Xq^ub obtained for the Auger 
data would have been obtained in the various models. Since there is a systematic uncertainty in the experimentally 
determined energy of the events (see fig. [T]) the table also lists the results obtained under the assumption that the 
energies of the events are systematically underestimated by 20%. 

Based on tabic IIIIl the distribution of the arrival directions of > 60 EcV cosmic-rays reported by the Auger 
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FIG. 8: The positions of the 27 Auger events with energy exceeding 5.7 x 10 eV, overlaid on the intensity map obtained in 
the biased model (see figure |4]). 





-^O , Auger 


Xi XuB Xb 


Pi Pub Pb 


Xc,ub{E > 57 EeV) 
Xc,ub{E > 68 EeV) 


2.7 
3.6 


0.7 1.8 
1.2 2.9 


0.9% (0.2%) 3.8% (2.1%) 20.4% (18.9%) 
1.5% (0.1%) 5.7% (2.7%) 27.6% (27.5%) 



TABLE III: Xc,ub obtained for the Auger data, denoted Xc, Auger, the average values of Xc,ub that are expected to be 
obtained in the various models (isotropic- Xj, unbiased- Xub, biased- Xb) for 27 events, and the probability that a value 
exceeding Xc. Auger would have been obtained in the various models (isotropic- Pj, unbiased- Pub, biased- Pb) for 27 events. 
The 2nd line gives the results under the assumption that the energies of the Auger events are systematically underestimated by 
20% (see fig. [TJ. Model parameters used: so = lO^'Mpc"^ (and so = 10"^Mpc"^ in brackets), s cx (l+z)^ , dlogn/dlog E ^ -2, 
and AE/E = 0.2. The galactic disk, |fo| < 12°, is excluded from the analysis. 

experiment is inconsistent with isotropy at ~ 98 % confidence level, and is consistent with a UHECR source distribution 
that traces the LSS. According to tablelTTl the probability of ruhng out isotropy at ~ 98% confidence with the observed 
number of events is not small. The observed distribution indicates a slight preference to a source distribution that is 
biased with respect to that of IRAS galaxies, that is, there is some indication that UHECR sources are more clustered 
than IRAS galaxies. As discussed in § |TT] and in § IIIIl the source number density, s{z = 0), may be constrained 
by the number of "repeaters", the number of sources producing multiple events. However, although a repeater may 
be present in the published Auger data, given the possible magnetic field deflections of few degrees (see § IIIBf) , the 
presence or absence of a repeater in the present data docs not provide stringent additional constraints on s(z = 0) 
(beyond sq > 10~^Mpc"^; see table |lT|. 

The following cautionary note should be made here. In their search for a correlation between UHECR arrival 
directions and the locations of V-C catalog AGNs, the authors of ref. [1^ have optimized several search parameters: 
the maximum angular separation of AGN and UHECR directions, the maximum distance of AGNs included, and the 
minimum energy of UHECRs included. While the angular separation and AGN distance parameters are irrelevant 
for our discussion, the minimum UHECR energy might be. As explained in the introduction, it is clear that the 
anisotropy signature should be searched for in the arrival distribution of UHECRs with energy threshold 5 x 10^^ eV. 
In ref. [l^ it was suggested to search for anisotropy with energy thresholds of 4 x 10^^ eV and 6 x 10^^ eV, and we 
have made a similar suggestion in sec. IIIIl based on table [Hi The arrival direction data that was made available 
by the Auger experiment has an energy threshold of « 6 x 10^^ eV, perfectly consistent with the a-priory energy 
thresholds suggested. However, the optimization of the energy threshold, to 5.7 x 10^^ eV, raises the concern that this 
optimization minimizes Pi, and hence leads to an overestimate of the significance with which isotropy is ruled out. 
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It is impossible to estimate the magnitude of such an effect based on the data that was made available. In order to 
confirm the conclusions drawn in the previous paragraph, one may search for the anisotropy signature in the Auger 
data using the lower energy threshold suggested, 4 x 10^^ eV (the probability of detecting a signal at this energy is 
larger than at higher energy, see table [n|, or repeat the analysis with an energy threshold of 5.7 x 10^^ eV using a 
larger data set (which will be accumulated on 1 yr time scale). 

V. SUMMARY 

We have derived the expected angular dependence of the UHECR intensity, under the assumption that the UHE- 
CRs are protons produced by extra-Galactic sources that trace the large scale distribution of luminous matter. 
All sky maps of the UHECR intensity, I{E, fl) with several energy thresholds E, are presented in figure [3] for 
an unbiased UHECR source distribution, where the UHECR source density is proportional to the LSS density, 
and in fig. [4] for a biased model, where the UHECR source density is proportional to the LSS density in over- 
dense region and vanishes elsewhere. A numerical representation of the intensity maps may be downloaded from 
http://www.weizmann.ac.il/~waxmaji/criso. The angular structure of the UHECR intensity map reflects the lo- 
cal large scale structure of the galaxy distribution, as can be seen by comparing figures [3] and |4] to figure [5l which 
presents the integrated galaxy density out to a distance of 75 Mpc. The anisotropy is larger for higher energy thresh- 
old, since the propagation distance increases at lower energy (see fig. [2]) and the contribution to the UHECR flux of 
sources beyond ~ 100 Mpc constitutes a roughly isotropic "background". We have used simple analytic arguments 
(§ IIIB[) to show that inter-galactic magnetic fields may modify the intensity map on scales of a few degrees, but not 
on larger scales (a result which is consistent with detailed semi-analytic and numeric analyses [3, E^)- 

We have defined a statistic, Xc (eq.[Tl]), that measures the correlation between the predicted and observed UHECR 
arrival direction distributions, and showed that it is more sensitive to the expected anisotropy signature than the power 
spectrum and the two point correlation function (see table [l|. The value of Xc for a given data set of UHECR arrival 
directions, the average value of Xq over realizations of the various models (isotropic, I, un-biased, UB, biased, B), 
and the distributions of Xq values expected in realizations of the various models in the limit of infinite UHECR 
source density, can all be straightforwardly calculated using the numerical representations of the UHECR maps at 
http://www.weizmann.ac.il/~waxman/criso. In order to avoid sensitivity to possible distortions of the UHECR 
intensity map by magnetic fields, we have used 6° x 6° bins in calculating Xc and excluded the Galactic plane region, 
\b\ < 12°. As can be seen by comparing results in tables [11 and HIl the exclusion of the Galactic plane region does not 
affect the results significantly. Table [Til demonstrates that the anisotropy signal is stronger at lower energy: Although 
the contrast of the fluctuations in the UHECR intensity is higher at high energy (see fig.[3l), the signal becomes weaker 
at higher energies since the number of observed UHECR drops rapidly with energy. At i? ~ 20 EeV a significant 
contamination from Galactic sources may be present and the deflection of cosmic-ray particles may become significant. 
In order to avoid distortions of the UHECR intensity map by these effects, it is advisable to choose an energy threshold 
of > 40 EeV. 

We have shown that the distribution of Xc is not sensitive to the assumed redshift evolution of the source density 
distribution, to variations in the intrinsic spectrum of protons produced by the sources, and to statistical errors in 
the experimental determination of event energies (see fig. [7] and table ITT)) . The distribution of Xc does depend on 
the assumed local (z = 0) source number density, sq: The average value of Xc obtained in realizations of the various 
models (I, UB, B) is independent of sq, but the width of the distribution is wider, and hence discrimination between 
models is more diffcult, for lower sq (see fig.[7J; and table HI)) . Thus, in order to determine the probability of obtaining 
a certain value of Xc in a given model (I, UB, B), one must make an assumption regarding the value of sq- As 
mentioned in § [Hi and discussed in detail in [l3], s{z = 0) may be constrained by the number of "repeaters", the 
number of sources producing multiple events. The last column of table HIl demonstrates that the presence or absence 
of repeaters may provide additional constraints on the source density (beyond sq > 10~^Mpc^'^) once the number of 
events observed above 40 EeV exceeds several hundreds. 

Figure [71i demonstrates that the distribution of Xc values, which is expected to be measured from the distribution 
of UHECR events with energies inferred by the experiment to exceed a certain threshold, is sensitive to uncertainties 
in the experimental calibration of the absolute energy scale: Since the anisotropy is stronger at higher energy, a 
systematic underestimate (overestimate) of event energies leads to a shift of the Xc distribution towards larger 
(smaller) values. Since the anisotropy is also stronger (at a given energy) for stronger bias, this implies that an 
underestimate (overestimate) of the absolute energy scale will lead to an overestimate (underestimate) of the bias of 
the UHECR source distribution. In order to distinguish between different bias models, the systematic uncertainty in 
the absolute energy calibration of the experiments should be reduced to well below the current ~ 25%. 

We have shown, using the Xc statistic, that the recently published > 5.7 x 10"'^^ eV Auger data are consistent 
with a source distribution that traces LSS, with some preference to an UHECR source distribution that is biased 
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with respect to the galaxy distribution, and inconsistent with isotropy at ~ 98% CL (see table HIH ). We have noted, 
however, that the optimization of the energy threshold made by the Auger collaboration analysis raises the concern 
that the significance with which isotropy is ruled out may be overestimated (see the end of § IIVP . In order to confirm 
our detection of a correlation with LSS, one may search for the anisotropy signature in the Auger data using a lower 
energy threshold, 4 x 10^^ eV, or repeat the analysis with an energy threshold of 5.7 x 10^^ eV using a larger data 
set (which will be accumulated on a 1 yr time scale). 

According to tabic |TT1 the detection of ~ 300 events above 4 x 10^^ cV is required in order to enable one to identify 
the predicted anisotropy signature, and to discriminate between different bias models, with statistical significance 
exceeding 99% CL. The experimental exposure required to accumulate this number of events is uncertain, due to 
the uncertainty in the absolute energy calibration of the UHECR experiments, which implies an uncertainty in the 
absolute normalization of the UHECR flux ([l^l, see also fig.[l|). Adopting, for example, the absolute energy scale of 
the HiRes experiment, the average number of events detected above 40 EeV is 300 per year for the Auger exposure 
of 7000 km^ sr, and about 50 events per year for the planned Telescope Array exposure of 1200 km^ sr [4^. These 
numbers are reduced by a factor of approximately 2 if the preliminary absolute energy scale of the Auger experiment 
is chosen (see fig. [ij . The numbers given in table |ll] imply also that the exposure accumulated within a few years of 
observation is unlikely to increase the significance of the detection to > 99.9% CL, unless the UHECR source density 
is comparable (or larger) than that of galaxies. 
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